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The invention relates to a magnetic resonance (MR) method for forming an 
.mage of an object from a plurality of signals sampled in a restricted homogeneity region of a 
mam magnet field of the magnetic resonance imaging apparatus, whereas a patient disposed 
on a table is moved continuously through the bore of the main magnet, such that an image is 
fonned over a region exceeding largely the restricted region of a main magnet field of the 
magnetic resonance imaging apparatus. 

The invention also relates to an MR apparatus and a computer program 
product for carrying out such a method. 

Since it is quite uncomfortable for a patient being disposed in a large bore of 
the main superconductive magnet of an MR imaging apparatus, there is a trend to use shorter 
magnets while moving the patient table or couch through the bore of the magnet. Such an 
apparatus is described e.g. in EP-A-1 024 371. m an MR imaging apparatus excitation pulses 
are applied to a restricted region of the magnet bore, in which the field is uniform. The data 
samples collected are Fourier transformed to form a volumetric image of the restricted 
15 region. A motor continuously moves a patient couch so that a region of interest passes 

through the region of good field. The collected data samples are corrected to compensate for 
the motion so that a volumetric image is formed of greater length than that of the restricted 
region. 

The MR imaging as mentioned above takes quite a long time since it is 
performed with conventional imaging sequences and even an overlap between successively 
imaged volumes is preferred in order to discard artefacts produced by errors in the central 
reg,on. Obviously, the r.f. transmit coil and the r.f. reC eive coil are disposed in the reference 
frame of the main magnet, i.e. in its direct vicinity. However, imaging by body coils disposed 
on different places of a patient was not an issue in this reference. Anyway, it is known that 
using body coils may cause additional artefacts which cannot be suppressed in an easy 
manner. 

Therefore, it is an object of the present invention to improve the MR imaging 
method with a relative short main magnet and a moving bed in such a manner, that imaging 
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will be performed in a much shorter time period. Another object of the present invention is to 
provuie an MR apparatus and a computer program product for carrying out the method. 

These objects are achieved by a method as claimed in Claim 1, by an MR 
apparatus as claimed in Claim 1 1 and by a computer program product as claimed in Claim 
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The invention has the main advantage, that shorter main magnets with a 
restricted homogeneity region can be used whereas data can be sampled in a much shorter 
time due to undersampling and therefore the scan time is largely reduced in relation to 
conventional MR imaging with a moving table or bed. Another advantage of the present 

10 invention is that the direction of the field-of-view can be oriented in any direction, which 
means that e.g. flow artefacts can be oriented in an appropriate manner or to orient single 
slices to some preferred orientation with respect to the tissue structure of the patient. 

These and other advantages of the invention are disclosed in the dependent 
claims and in the following description in which an exemplified embodiment of the invention 

15 is described with respect to the accompanying drawings. Therein shows: 
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Fig. 1 diagrammatically a magnetic resonance imaging system in which the 
present invention is used, 

Fig. 2 an example of a acquisition scheme, 

Fig. 3 an acquisition scheme in k-space with a moving bed, and 

Fig. 4 a simplified presentation for explanation of the present invention. 

Basic pri nciples and notations 

In Fig. 1 a patient 1 is placed fixedly to a moving table or bed 2. A set of local 
coils LC, to LC 5 is located fixedly to the patient 1, thus the arrangement "patient 1" "bed 2" 
and "local coils LC, to LC 5 " are disposed in a fixed relationship to each other in a coordinate 
system 3 which is related to the patient. A vector x denotes the position relative to this 
coordinate system or frame of reference 3. Further a main magnet 4 is provided with a bore 5, 
in which the bed 2 can be moved to and fro in direction of arrow 6. The magnet/gradient 
system 7 comprises the main magnet 4 and a number of receive coils 8 - generally one or 
none - and a transmitter coil (not depicted). These receive coils 8 and transmit coils are 
denoted as being "global". The elements of the magnet/gradient system 7 are all related in 
another fixed relationship which is defined by another coordinate system 9 related to the 
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main magnet 4. It is further assumed that the receive coils 8 and the transmit coils are large 
and have smooth sensitivity patterns in the direction of the bed motion (arrow 6). The magnet 
coordinate system 9 is offset by a distance A relative to the patient coordinate system 3. The 
property of the patient 1, e.g. the proton density of his body, is obviously given in the patient 
coordmate system 3, expressed as pg) . Ideally, the same holds for the receive-sensitivity 
pattern S, (*) of one of the local coils LQ, where i = 1 to 5. On the other hand, the sensitivity 
pattern of the receive coil 7 is written as g/(x - A) . 

Now the expression "receive situation" is introduced, which means a block of 
measurements contiguous in time, during which the table 2 has travelled a relatively small 
distance (e.g. one half or one third or something in between) of the "homogeneity size" of the 
mam magnet 4. During a single receive situation all magnetisation preparations (like 
presaturation slabs, volume selections etc.) "travel" in principle with the patient 1 Also 
durmg a single receive situation, all of the interesting k-space is covered, except for any 
(potentially significant) undersampling. The full scan over k-space covers K receive 
stations. In principle, the MR acquisition sequence is identical over all receive situations 
except that the "magnetisation preparations" may change place and that there may be some' 
offset in the acquired k-space grid. This is however a refinement of the concept. 

During any receive situation, every coil LQ to LQ receives information from 
the parent 1. And so (I + J) sets of information are delivered, where I is the number of local 
recerve coils and J is the number of global receive coils. Some of these sets may be 
practically zero, e.g. if a receive coil LC, to LC 5 is completely outside the homogeneity 
volume of main magnet 4. There are K receive situation, Thus, in total, information has been 
gathered m K * (1 + J) "receive instances". 

Basic idea of the invgntirm 

The clinical doctor wants information from a relatively large volume of 
interest of the patient 1. Data is acquired while the table or bed 2 is moving relative to the 
main magnet 4. The displacement A thereof relative to the magnet 4 is known at any time 
The whole of the MR acquisition is split up in a number of receive situations. During each 
Situation, the excitation and presaturation profiles "travel" with the patient 1 . The acquired 
data are phase-corrected for the offset of table 2 at which they are acquired. During a receive 
situation, information from the patient 1 is acquired using any MR imaging sequence of any 
onentation. In practice it may be useful either to orient artefacts (e.g. flow artefacts) in an 
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appropriate direction, or to orient single slices to some preferred orientation with respect to a 
specific human-tissue structure. However, the information may be seriously undersampled. 
An example of this method is sketched in Fig. 2. For every receive situation, we acquire, 
nominally the volume indicated by a box 1 1. This means that a multitude of points 12 will 
fold onto the results of this receive situation. This multitude may be a one-dimensional row, a 
two-dimensional row or even a continuum of folding locations which is the case in e.g. spiral 
acquisition. 

Now, every receive situation has the same subsampling pattern. The only 
difference between different receive situations is the weighing with respect to the properties 
of following selected factors: 

1. the excitation or saturation profile (with respect to the patient 1), 

2. the transmit pattern of the transmit coil, 

3. the receive patterns of the fixed receive coils 8, and 

4. the frequency-response pattern of the receiver system. 

Concerning items 2 and 3 it is remarked that in theory, with a moving table 2, 
these patterns varies also within a situation - but if the motion per situation is not exceedingly 
large, this may be neglected. 

In addition, every receive coil (whether local or global) has its own sensitivity 
pattern. So every receive instance, which is the combination of a receive coil and a receive 
situation, "sees" a different "overall pattern", which is a product of the coil-sensitivity pattern 
and the combined pattern (excitation, saturation, transmit and frequency-response) of the 
receive situation. 

The multitude of folding points 12 is unfolded using all the "overall patterns" 
of every receive instance. This is done in a manner as in the SENSE method (see e.g. K. 
Pruessmann et. al. in Proc. ISMRM, 1998, abstracts pp. 579, 799, 803 and 2087) - although, 
formally, the unfolding matrix may be larger. Nevertheless, the system of equations may be 
reasonably stable: e.g. for each situation, many folding points 12 may practically "see" a zero 
pattern, because they fall outside the excitation slab. It is further to be noted that the table 
motion is flexible: it may be linear, non-linear and even non-monotonic or two-dimensional. 
The restriction is that it must not travel excessively during one receive situation, e.g. no more 
as half the homogeneity volume. 
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Measuring the patterns 
Calibration 

In SENSE, a map of coil sensitivity profiles must be available. Such a map can 
be provided in a reference or coarse calibration measurement with a well-defined phantom. 
Another possibility is using a low-resolution map of the ratio of signals from an array coil 
element and from a reference coil of the patient on the table. In the present invention a 
similar reference measurement is provided. Although data can generally be acquired at any 
table-motion pattern, including two-dimensional table motion, here a linear one-dimensional 
movement will be used for data acquisition. 

For a better understanding of the involved principles, first a review: 

1. At each position of the patient (x ) some position-dependent tissue-properties 
(all data about proton density, T u T 2 , flow, diffusion...) are known 

2. Also at each position of the patient ( x ) the sensitivity pattern S, (x) of the 

local coils LCi can be obtained. Although there might be some dependence on A due to 
coupling to global coils, this is neglected here. 

3 • Tne pattern of all the coils (global or local, transmit or receive) may be 

influenced by dielectric resonances, particularly at main magnet fields of 3T and higher. It is 
assumed that this effect can be seen as purely multiplicative on the coil sensitivity of all coils 
- an effect that is only dependent on the patient and the position therein, so not on the coil or 
coil-position. For a main magnet field of 1,5T the effect is negligible. The dielectric 
resonance effect is denoted as D(x) . Each coil actually "sees" with a sensitivity D(x) • S,(x) . 
4. The global coils also have their sensitivity pattern, which is approximated as 

being independent on the actual position of the patient within the magnet frame (except, of 
course, for dielectric resonance effect, which is explicitly accounted for). This sensitivity 
pattern is expressed as (x + A) . 

5 - 1116 total receive-pattern also may contain effects of receiver frequency 

response, denoted as F(x + A) at a given readout gradient direction. 

For a given measurement type, the transverse magnetisation is a non-linear 
effect depending on both the transmit pattern - Le. on the position in the magnet coordinate 
system 9, if a global coil is used - and the tissue property - i.e. on the position in the patient 
coordinate system 3. Since this is very difficult to handle, an approximation is necessary, 
which can be achieved by using an FFE measurement with very small tip angles (about 5 
degrees or less). In that case, the transverse magnetisation can be approximated as 
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/(*).«(* + A),where /(*) is the "spin density" function, which is a tissue property for a 
given measurement (in the proposed case, it mainly depends on proton density), a is the tip 
angle, which is approximated as being proportional to D(x) • g lramtnlt • (3c + A) . 

5 Problems 

For the reconstruction of the proposed method as described in more detail 
below, ideally one would like to know all information of *,<*), gj( x + A ) and D(Jc) (but 
the last one is not crucial). The value of F(x + A) is assumed to be known beforehand, and 
the value of /(*) is (at least in principle) not an interesting result of the calibration 
10 measurement. 

The problem is that not all these parameters are measured. For a given table 
position A, the following is measured: 

- the signal from local coils: F(x + A)-S, (xyDffi.g^^ (x + A).D(x).f (x) . 

- the signal from global coils: F(x + A)-gj(x + A).D(x).g lnmmll (x + A).D(x).f (x) . 

A plain division of those two signals results in S t (x)/ gj(x + A) - which is of no direct use if 
gj(x + A) is not uniform. 

First solution 

The patterns of the global coils are assumed to be known ( i.e. at least one of 
them). In practice, that function may be load-dependent, but this can be overcome by the fact 
that there is a calibrated load-dependency of the sensitivity function, and the load can be 
easily measured. 

A number of discrete reference-scan segments are acquired. During such a 
segment, the table stands still. Each segment results in a full "image". The table is then 
stepped over a number of segments, so for a number of different values of A . The 
knowledge of gj(x + A) allows calculating the matrix S, (3c). Mathematically, each segment 
would give the result over the full extent of x , but in practice, each segment will provide 
accurate results only over a sub-range of all x . These partially overlapping ranges can be 
combined using least-squares fits, resulting in the full map of S,(x) . 
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Second solution 

The three-dimensional measurement is oriented in such a way that the 
frequency-encoding is parallel to the table motion direction - the z-direction. The sampling 
bandwidth is very high in that direction. The other two directions are phase-encodings, one of 
winch is "fast" compared to table motion. The "slow" phase-encoding direction is called k y 
In the patient coordinate system 9 the main magnet 4 is slowly moving with respect to the 
patient During the motion of the magnet, profiles are acquired with linearly increasing ky In 
Fig. 3 a line segment 13 indicates the portion of the patient 1 that is detected by the 
homogeneity volume of the main magnet 4. 

L is the distance that the main magnet 4 covers for a full set of phase-encoding 
profiles. The table speed is arranged in such a way that L is maximally half of the 
homogeneity volume - but less (e.g., 1/4 ) is preferable. In Fig. 3 L is chosen at 1/2 The 
consequence is that every profile has been sampled for at least two different positions of the 
mam magnet 4 (except for the outer edge regions, which are disregarded here). For example a 
profile k y has been measured for both the displacement A, and for the displacement A 2 This 
allows to approximate the profile that would have been measured if the main magnet 4 were 
at position A n/ (as long as this value is in between A, and A 2 ). This can be done by 
interpolation (e.g. a linear interpolation if two profiles are available, or a cubic spline if there 
are four profiles). That can be done for every profile, which allows a complete reconstruction 
of the sensitivity-calibration measurement at magnet-offset A re/ . 

However, this can also be done for any value of A*/ . That principle allows 

acquiring: 

- the signal from local coils: F ( x + A)- S> (*) . m . g _ (5 + A) . ^ . /( - } 
-the signal from global coils: F(x + A). gj(x+ A W xy g ^ 

This holds, in contrary to the "problem section, not only for one value of A but for all values 
thereof. This allows to: 

- separate the signal from local coils into F(x + A) • g lnmmll (x + A) and 

- separate the signal from global coils into F(x + A)-gj(x + A). gtnmmll (x + A) and 
D(x).D(x).f(x). 
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This allows the derivation of S,x and of gj(X + A) . Actually, if F is known, then also 

**— + A) could be derived, but unfortunately not the dielectric resonance D ( 3c) - but 
that is not crucial for the sequel. 

Receiving and reconstructing th* art,.*) A m 

In the section "basic idea of the invention" the receiving and reconstruction of 
data has been explained to a great extent. The data is acquired using relatively large steps 
between profiles, resulting in a (nominally) small "folding volume", i.e., into lots of folding 
For non-cartesian sequences this means that data is sampled relatively sparse. The full extent 
of k-space is acquired during one scan situation, and the displacement during one scan 
situation is a fraction (e.g. 1/3) of the homogeneity volume. Every incoming sample is 
acquired at a very specific value of A,^ . This is taken into account by multiplying every 
sample of the incoming data with exp(- ikTA sample ) . Here, T is the coordination 
transformation from patient coordinates to scan coordinates, which is relevant for oblique 
scanning. That operation "displaces" the acquired data to the centre of the patient, even if that 
is way out of the homogeneity volume. Obviously, one can simplify things by correcting with 
a fixed Apro/ue per profile. 

The scan situation has been acquired at an average offset A w . This is relevant 

to estimate: 

- the transmit profile for that situation: g\ ,(x) = e fS + X ^ 

« transmit \*J — & transmit K x + &ave ) > 

- similarly for the global-coils receive profile, and 

- similarly for the frequency response function for that situation. 

The excitation profiles and presaturation profiles are fixed to the patient 
coordinate system, i.e. travelling with the table 2 as seen from the main magnet 4. So these 
are known in x . Also the local coil profiles are known in x (see calibration section). All 
these coil profiles can be multiplied into an overall receive pattern for each coil. In total, if 
there are K receive situations, then information has been gathered in K°(I + J) receive 
instances. Each instance has a different overall receive pattern. This allows for a SENSE-like 
reconstruction, if there are less than K°(I + J) folded points over the entire patient 1. In 
principle, the acquisition stays the same over all receive situations. As a refinement, the 
phase encoding may be offset by a small amount. This will give some extra phase-encoding, 
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on top of all other mentioned encodings, which will improve the stability of the SENSE- 
reconstruction. 

Geomet ric correction 

All of above reasoning assumes a perfectly linear gradient system. It assumes 
that a table-displacement of A will cause (a) some different weighting, which has been 
extensively treated above, and (b) a phase modification of exp(/^A_ /e ) to a sample. 
Unfortunately, geometric distortion is another significant complicating factor. An image 
acquired at offset A, is geometrically distorted when compared to an image acquired at offset 
A 2 , even if all weightings and the displacement ( A 2 -A, ) have been accounted for. 

This problem can be solved since the distortions tend to be small where the 
excitation and transmit profiles tend to be large. This can be accomplished (a) by making 
sequences in such a way that there is not too much excitation in geometrically 
inhomogeneous areas, and (b) by designing the system with low transmit-sensitivity in these 
areas. The solution is an iterative reconstruction, the explanation thereof is given visually in 
Fig. 4. The main idea is, in essence, to first reconstruct the resulting image by disregarding 
the presence of geometric distortion. Of this result, all image regions are distorted in as many 
ways as there are scan situations. The distorted minus the undistorted versions provide an 
estimate of "folding distortion-errors". If the contributions (or weights) of every receive- 
instance to the final image are known, the estimate of the folding distortion-errors can be 
extracted from the result. 

In Fig. 4 (a) a very simplified patient 15 is depicted which diagrammatically 
consists of two parts at different positions: a first part shown as block 16 and a second part 
shown as heart 17. Fig. 4 (b) shows the imaging while the block 16 is largely in the centre of 
the main magnet and the heart 18 is more remote and thus distorted. As shown in Fig 4 (c) 
the distorted heart 18 folds onto the block 16, due to undersampling. Yet, the heart 18 is 
fainter, since it has seen less excitation. In another receive-situation according to Fig 4 (d) 
the heart 17 is in the centre ofthe main magnet 4 and the block 19 is more remote and thus 
distorted. Due to undersampling, this is measured in that receive-situation, as shown in Fig 4 
(e). So the images in Fig. 4 (c) and Fig. 4 (e) are actually measured, Fig. 4 (f) shows the 
result after SENSE-unfolding, which consists of two linear combinations of Fig. 4 (c) and 
Fig. 4 (e). It is imperfect, e.g. the right part of it consists ofthe block 16 plus an unsubtracted 
heart 18 plus a "wrongly subtracted" heart 17 (dashed) and a second-order error (dotted). In 
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Fig. 4 (g) is shown a purposely distortion in the same manner as in (b), under the condition 
that (t) is available. Fig. 4 (h) shows the final result where part of (f) is subtracted from its 
folded counterpart, and part of (g) is added, resulting in only second-order artefacts. 



